Electron-electron interactions in decoupled graphene layers 
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Multi-layer graphene on the carbon face of silicon carbide is an intriguing electronic system which 
typically consists of a stack of ten or more layers. Rotational stacking faults in this system dramati- 
cally reduce inter-layer coherence. In this article we report on the influence of inter-layer interactions, 
which remain strong even when coherence is negligible, on the Fermi liquid properties of charged 
graphene layers. We find that inter-layer interactions increase the magnitudes of correlation energies 
and decrease quasiparticle velocities, even when remote-layer carrier densities are small, and that 
they lessen the influence of exchange and correlation on the distribution of carriers across layers. 
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I. INTRODUCTION 

Graphene layers prepared by the thermal decomposi- 
tion of silicon carbide (SiC^PEl might play a pivotal role 
in carbon-based analog electronics because of their high 
carrier mobilities, and because of the possibility of direct 
growth on a well understood single-crystal semiconductor 
substrate. This form of graphene, usually termed epitax- 
ial graphene, can be grown directly on either silicon or 
carbon terminated faces of SiC and is scalable to large 
area circuits. 

Epitaxial graphene on the carbon face typically con- 
sists of tens of graphene layers (see however Ref. [T0|h 
with rotational stacking faults that lead to an electronic 
structure which is practically indistinguishable^^ from 
that of a sample with many electrically isolated single- 
layer graphene (SLG)P^ layers. It is commonly be- 
li eve( pmili31i6J t h at t h e i ayer c i sest to the SiC typi- 
cally has a rather high carrier density « 10 12 cm -2 , due 
to charge transfer from the substrate, and that the over- 
layers are nearly neutral and unimportant for most ob- 
servables. Even though multi-layer graphene (MLG) on 
the C-face of SiC has been shown to possess many of the 
distinctive features ofSLG, including the half-quantized 
quantum Hall effectP^ 2 ^, two key questions arise: i) to 
what extent do Coulomb interactions between electrons 
in nearby layers distinguish the electronic properties of 
MLG layers from those of SLG? ii) what is the role 
of inter-layer Coulomb interactions in determining the 
electron-density distribution which achieves equilibrium 
between layers? These are the main issues we address in 
this article. 

The system we study is sketched in Fig. [I] We as- 
sume that the graphene layers are perfectly decoupled 
from the point of view of single-particle tunneling, but 
take Coulomb interactions between electrons in differ- 
ent layers, which provide a source of (two-particle) cou- 
pling, fully into account. We study, (i) the impact of 
intra- and inter-layer Coulomb interactions on the Fermi- 
liquid properties of MLG two-dimensional electron sys- 
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FIG. 1: (Color online) A stack of £ = 1 ... N graphene layers 
is placed between top and bottom dielectrics. The layers are 
coupled solely by inter-layer Coulomb interactions (red wiggly 
lines). The dielectric constant in the region above the ^-th 
layer is labeled by ei , while the one of the substrate is labeled 
by ejv+1- The separation between the £ + 1-th and the £-ih 
layer is labeled by c^+i^. 



terns, and (ii) the combined role of Hartree, exchange, 
and correlation potentials in determining the equilibrium 
distribution of carrier densities across layers when charge 
is transferred to the many-layer system from the SiC sub- 
strate. The theory sketched below applies to MLG on the 
C-face of SiC with an arbitrary number N of layers, to 
the decoupled layers sometimes found on the surface of 
bulk graphite^, and to the experiments by Schmidt et 
alW* on micromechanically exfoliated decoupled layers. 
For the sake of definiteness and simplicity, we present 
explicit numerical results only for n-type carriers and for 
the N = 2 (double- layer) case. 

This manuscript is organized as follows. In Sect. [TT] 
we present the model we have used to describe tunnel- 
decoupled graphene layers in the presence of intra- and 
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inter-layer electron-electron interactions, and discuss the 
linear-response functions which control the physical prop- 
erties in which we are interested. In Sect. lIIII we describe 
the minimal microscopic theory that allows us to calcu- 
late quasiparticle velocities, ground-state energies, and 
chemical potentials. In Sect. IV we present a general 



scheme to tackle the electro-chemical equilibrium prob- 
lem in MLG. Finally, in Sect.[V|we present our main nu- 
merical results for double-layer graphene. Sect. |Vl| con- 
tains a summary of our main conclusions. Appendix [A] 
collects some lengthy expressions for the exchange and 
correlation energies of double-layer graphene that are im- 
portant for the technical details of our calculations. 



II. MODEL HAMILTONIAN AND 
LINEAR-RESPONSE THEORY 

Our model Hamiltonian contains massless-Dirac- 
fermion kinetic-energy terms and intra- and inter-layer 
Coulomb interactions: 
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Here v is the bare Fermi velocity, taken to be the same 
in all the £ — 1 ... N tunnel-decoupled layers, S is the 
area of each layer, Vw (o) is the matrix of bare Coulomb 
potentials, and 
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is the density-operator for the £-th layer. Greek let- 
ters are honeycomb-sublattice-pseudospin labels and a = 
(a x ,cr y ) is a vector of Pauli matrices. We also in- 



troduce the coupling constant!^ 



- .2 



e 2 /(hv), whose 



value is « 2.2 if for v we use the SLG Fermi velocity 
vp ~ 10 6 m/s. 

Several important many-body properties of the Hamil- 
tonian H are completely determined by the N x TV sym- 
metric matrix x(g,o;) whose elements are the density- 
density linear-response functions 
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with ((AjB})^ the usual Kubo productP. Within the 
random phase approximation (RPA) these functions sat- 
isfy the following matrix equation, 



x 1 (Qiu) = x 1 (q^) ~ V(q) 



(4) 



where Xo(q^) is a N x N diagonal matrix whose ele- 
ments a;) are the well-knowrP3MU noninteracting 
(Lindhard) response functions of each graphene layer at 
arbitrary doping ri£ = N^/S. The off-diagonal (diago- 
nal) elements of the matrix V = {V^'}^'=i...jv represent 
inter-layer (intra-layer) Coulomb interactions. 



III. QUASIPARTICLE VELOCITIES, 
GROUND-STATE ENERGY, AND CHEMICAL 
POTENTIALS 

The Fermi-liquid parameters of MLG can be calcu- 
lated from the knowledge of the retarded quasiparticle 
self-energy E^. Our result s are based on the so-called 
"GoW" approximatio n 23 * 27 " 30 * in which the self-energy T,£ 
is expanded to first order in the dynamically screened ef- 
fective interaction W^H 



W(q,w) = V(q) + V(q) X (q,oj)V(q) 
= [V- 1 (q)-Xo(q,^]- 1 ■ 



(5) 



In this paper we limit our attention to many-body quan- 
tities that can be expressed solely in terms of integrals 
along the imaginary frequency axis where the Lindhard 
function has a smooth dependence on its parameters. 

The microscopic expression for Tig in terms of W can 
be obtained by a straightforward generalization of the 
theory of Ref. [29] to a multicomponent system (h = 1): 
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where f3 (k B T) 1 , k' = k + q and uo' n m = uo n + Q rn . 
In Eq. ^ uj n = (2n + l)n/P is a fermionic Matsub- 
ara frequency while the sum runs over all the bosonic 
Matsubara frequencies Q m = 2m7r//3. The factor in 
square brackets in Eq. (|6|, which depends on the angle 
6k,k+q between k and k + q, captures the dependence of 
Coulomb scattering on the relative chirality s of the in- 
teracting electrons. The Green's function G^l(k,iu) = 
l/[iuj — describes the free propogation of states 

with wavevector fc, Dirac energy & }S (fc) = svk — \i£ (rel- 
ative to the chemical potential fi£ of the £-th layer) and 
chirality s = =b. Note that in the r.h.s. of Eq. (|6| there 

are no terms involving products of the form WwG^\ 
with £' ^ £. The reason is that, due to the absence 
of inter-layer tunneling, bare propagators are diagonal 
in the layer index: in the diagrammatic language (see 
Fig. 8.17 in Ref. l23|) . screened interaction Wu' wavy lines 
closed on a bare propagator cannot begin on layer label 
£ and terminate on layer label £' ■=£ £. 

After analytic continuation from imaginary to real fre- 
quencies, iuj — » uj + ir], the renormalized Fermi velocity v\ 
for the quasiparticles in the £-Xhi layer can be expressed in 
terms of the wavevector and frequency derivatives of the 
retarted self-energy S^ et (/c, uo) evaluated at the Fermi sur- 
face (k = fcp^) an d a ^ the quasiparticle pole uo = £g j+ (fc): 



(7) 



l-^e£f(fc, W )| fc=fcF ^ =0 
As explained elsewher e 1 23 * 28 * 29 ^, the derivatives in Eq. ^ 
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can both be expressed in terms of integrals along the 
imaginary frequency axis. 

The ground-state energy of MLG can be easily evalu- 
ated from T~L. Apart from the trivial kinetic energy con- 
tribution fekin related to the first term in T-L, it contains 
intra- and inter-layer interaction energy contributions. 
To evaluate th e inte raction energy we have followed a fa- 
miliar strateg y 23 ! 24 * by combining the coupling-constant 
integration algorithm with the fluctuation-dissipation 
theorenPEI] Following Ref. HQ we choose the total en- 
ergy of undoped MLG (Fermi energy at the neutrality 
point in all layers) as our zero of energy. We separate 
the contribution that is first order in a ee , the exchange 
energy, from the higher order contributions convention- 
ally referred to in electron-gas theory as the correlation 
energy. Generalizing the theory for SLGP^I to a multi- 
component system, we end up with rather cumbersome 
expressions for the exchange fe x and RPA correlation en- 
ergies Ss^ PA of MLG (per electron) measured from the 
reference undoped system. Explicit expressions for the 
N = 2 case have been reported in Appendix |Aj Because 
these expressions involve only imaginary axis frequency 
integrals they are easy to evaluate accurately by quadra- 
ture. 

The determination of the electro-chemical equilibrium, 
discussed in more detail below, requires knowledge of ex- 
change and correlations contributions to the chemical po- 
tentials fi£. These can be easily calculated from 

d(nSe tot ) 

where n = ni+n2 + . . . is the total density and fc t ot = 
fekin + fe x + Se^ PA is the total energy per electron, a 
function of 774, ... , n^. 

IV. ELECTRO-CHEMICAL EQUILIBRIUM 

The equilibrium density distribution across a multi- 
layer is achieved when the electro-chemical potentials in 
each distinct electronic region are identical. For a TV- 
layer graphene system grown on a (SiC, say) substrate we 
identify iV + 1 electronic regions; the additional system is 
a buffer layer (not shown in Fig. [I]) between the bulk sub- 
strate and the graphene layers which is positively charged 
as a consequence of charge transfer to the graphene sys- 
tem. The equilibrium densities in all the N + 1 layers 
are determined by N + 1 equations which express the 
discontinuity of the electric dispacement across charged 
layers. With the notation introduced in Fig. [I] Gauss's 
law implies that for £ = TV, N — 1, . . . , 1 (from bottom to 
top) 

ei+iEi+i - eiEi = -47rea^ , (9) 

where <jg is the areal electron density in the £-th layer. We 
can view E\ as a quantity which is determined by a gate 
voltage. As long as the distance dxB from the top gate 



electrode to the MLG system is large, "chemical" poten- 
tial contributions are negligible so that eEid^G = Vtg- 
We can write down an equation similar to ([9| for the 
buffer layer, 6at+2^at+2 — £n+iEn+i — +47re<Jb. Again 
we can view E^+2 as an experimentally controllable pa- 
rameter (fixed by a bottom gate, e^Ar+2^BG = Vbg)- 
Here <Jb is the density of positive charges in the buffer 
layer. This model assumes that the substrate is essen- 
tially insulating. In the absence of gates or unintended 
dopants that might create electric fields, we would nor- 
mally expect Ei = E^+2 = 0. In this case we obtain the 
charge neutrality condition: Y^z=i G t = °~h- In general 
<7b should not be considered a free parameter. It is de- 
termined by the total graphene charge and the difference 
between the two experimentally fixed electric fields Ei 
and En+2- For MLG samples grown on SiC depends 
on growth parameters in a way which is at present not 
well understood. 

As explained above, our convention for the zero-of- 
energy of chemical energy of each layer is that is is zero 
at the the Dirac point, i.e., at electrical neutrality. It is 
convenient to also choose an explicit global zero for the 
electric potential, which we take to be its value on the top 
layer. Given the charge densities, the electrical potentials 
in each layer can be conveniently calculated iteratively 
starting from the top layer. For £=l,2,...,iV — lwe 
have that 

V i+1 =V i + eE i+1 di+ li i . (10) 

Here di+i^ is the separation between the £ + 1-th and 
the £-th layer. The condition for equilibrium between a 
graphene layer and the layer below it is: 

V>t+l + Vg+l = fl£ + Vi . (11) 

We need one more equation to fix the densities and that 
is the equilibrium condition between the bottom layer 
and the buffer: + (j) = Vn + Mat- Here <j) represents the 
microscopic physics which causes electrons to spill out of 
the buffer layer (where they are likely poorly bonded). It 
is known that <p is sensitive to the arrangement of carbon 
atoms in the buffer layer. If all of the graphene layers 
are negatively charged there will be a large electric field 
between the N-th layer and the buffer layer whose sign 
will tend to repel electrons, i.e. to make Vb < Vn- The 
value of (j) must therefore be positive and it should be 
chosen to match experimental results for the total charge 
of all layers in a MLG systems. We assume that (j) is 
fixed once a sample has been prepared, i.e. that it is 
not influenced by gate voltages, external magnetic fields, 
or other parameters that are routinely used to alter the 
electrical properties of two dimensional electron systems. 

V. NUMERICAL RESULTS AND DISCUSSION 

We now turn to the presentation of detailed illustra- 
tive numerical results for double-layer graphene (DLG). 
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The bare intra- and inter-layer Coulomb interactions are 
influenced by the layered dielectric environment. For the 
N = 2 case, a routine electrostatics calculation^ implies 
that the Coulomb interaction in the £ — 1 (top) layer is 
given by 



Vn(q) 



47re 2 



[(e 2 + e 3 )e" d + (e 2 - e 3 )e-" d ] , (12) 



qd] 



where 



D(q) = [(e 1 +e 2 )(e 2 + e 3 )e« d + (e 1 -e 2 )(e 2 -e 3 )e-i d ] (13) 

and d is a shorthand notation for the inter-layer dis- 
tance c?2,i. The Coulomb interaction in the bottom layer, 
^22(^)5 can be simply obtained from V\\{q) by inter- 
changing 63 <o> e\. Finally, the inter-layer Coulomb in- 
teraction is given by 



V 12 (q) = V 21 (q) 



87re 2 
qD(q) 
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(14) 



In the case N = 2 the matrix equation Q can be easily 
inverted and the screened potentials W in Eq. ([5| can 
be written in a particularly compact form 3 - 2 ^ 



Wu(q,uj) 
and 



Vn(q) + \V? 2 (q) - V u (q)V 22 (q)}x 2 0) (q, w) 



W 12 (q) 



e(q,u) 
e(q,u) ' 



(15) 
(16) 



The expression for W22(q) can be obtained from Eq. (15) 
by interchanging 1 <H> 2. In Eqs. (15)-(16) we have intro- 
duced the dielectric function 



S{q,w) = [l-V 11 (q) X ( i\q,io)}[l-V 22 (q)x 2 l "(q,io)] 



(0), 



V? 2 (q)x { °\q,u)x 2 U '(q,u) . 



(0), 



(17) 



Exchange and correlation energies in MLG systems de- 
pend both on interactions on the Fermi wavelength scale 
which influence correlations between carriers and on in- 
teractions at shorter length scales which influence cor- 
relation with the Dirac sea background. This contrasts 
sharply with the case of ordinary two-dimensional elec- 
tron gases for which many-body properties depend only 
on interactions on the Fermi wavelength scale. Because 
MLG layers are separated by atomic length scales and 
carrier densities are always small on atomic scales /cp id 
is typically small. Taking the limit qd — »> in Eqs. ( 12 )- 



(14) we find that the typical carrier-carrier interactions 
are layer independent with effective dielectric constant 
(^i + e 3)/2- Similarly, taking the limit qd — » oo we find 
that carrier-background interactions are approximately 
independent in different layers, and that they have an 
effective dielectric constant determined mostly by imme- 
diately adjacent layers. 



With these explicit expressions in our hands we are in 
the position to calculate the quasiparticle self-energies 
and velocities v\, the exchange fe x and RPA correlation 
fe^ PA energies, and to solve the electro-chemical equi- 
librium problem outlined above. For the presentation of 
the numerical results for the double-layer system we in- 
troduce the density polarization £ = (n^ — n\)/n: ( = 1 
when the carrier density is non-zero only in the bottom 
layer, while ( = when the two layers have identical 
carrier densities. 

In Fig. [2] we present results for the quasiparticle ve- 
locity enhancements in the low and high density layers, 
v\jv and v\jv, as functions of their carrier density for 
several different values of the density in the opposite 
layer. The quasiparticle veloc ity i n a single-layer has a 
large enhancement due mainl y 29 ! 33 * to exchange interac- 
tions between the carriers and the Dirac sea. Carrier- 
carrier interactions yield positive contributions to both 
the wavevector and frequency dependence of the self- 
energy, which partially cancel as contributions to the 
renormalized quasiparticle velocity much as they do in 
the ordinary two-dimensional electron gas. The end re- 
sult is that the enhanced velocity due to exchange in- 
teractions with the Dirac sea largely survives, especially 
at low carrier densities. The parameters of this calcula- 
tion were chosen with DLG on SiC in mind. The general 
result is that enhanced quasiparticle velocity again sur- 
vives, with quantitative changes due to electron-electron 
interactions. Two aspects of these results may appear 
counter-intuitive at first glance. First of all, we see that 
inter-layer interactions are important even when the car- 
rier density in remote layers is zero. Because the Dirac 
bands are gapless and because both bands have 7r-orbital 
character, their wavevector and frequency dependent po- 
larization function is important in forming the dynami- 
cally screened Coulomb interaction even in the absence 
of carriers. Secondly, we notice that increasing the den- 
sity of a remote layer does not necessarily weaken the 
screened interaction in the layer of interest. The ori- 
gin of this property can be traced to the second term in 
Eq. (15) which is not intuitive and captures the mutual 



screening response of double-layers including inter-layer 
interaction effects. 

In Fig. [3] we present numerical results for the DLG 
exchange and RPA correlation energies per carrier, fe x 
and fe^ PA . In this case we plot energies as a func- 
tion of total carrier density for several different layer po- 
larizations in order to illustrate a relationship to well 
known dependences on spin-polarization that we di scuss 
below. Exchange energies are positive because^^ they 
are calculated relative to zero carrier density using the 
Dirac point self-energy of this limit as the zero of energy. 
The increase in exchange energy with carrier density in 
graphene has the physical consequence that corrections 
to the RPA are expected^ to enhance screening, instead 
of weakening it as in an ordinary two-dimensional elec- 
tron gas. Notice that the exchange energy of DLG is 
larger than that of SLG. The two are equal for £ = 1 
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FIG. 2: (Color online) Quasiparticle velocities in double-layer 
graphene on SiC [ei = €2 = 1.0 and 63 = 6.6 (SiC dielectric 
constant)]. The data shown in this figure have been obtained 
by setting d — 3.35 A and a ee = 2.2. Panel a) Quasiparticle 
velocity v* in the top layer (in units of the bare velocity v) as a 
function of the density m in that layer (in units of 10 12 cm -2 ), 
for different values of the density n 2 in the opposite layer. 
Panel b) Quasiparticle velocity v% in the bottom layer (in 
units of the bare velocity v) as a function of the density 712 
in that layer (in units of 10 12 cm -2 ), for different values of 
the density m in the opposite layer. Circles label data for the 
quasiparticle velocity in single-layer graphene^ on SiC. 
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FIG. 3: (Color online) Interaction energies in double-layer 
graphene on SiC, measured from the reference undoped sys- 
tem. All parameters are the same as in Fig. [2] Panel a) 
The exchange energy Se x in units of the Fermi energy sf as a 
function of n (in units of 10 12 cm -2 ) and for different values 
of C t . Panel b) Same as in the top panel but for the RPA 
correlation energy Se^ PA . Crosses label the RPA correlation 
energy for £ = 1 if the inter- layer interaction V12 is set to 
zero. Circles label the exchange and RPA correlation energies 
in single-layer graphene^ on SiC. Note that the correlation 
energy obtained neglecting inter-layer interactions (crosses) is 
practically equal to the single- layer graphene result (circles). 



because the exchange energy depends only on intra-layer 
interactions [see Eq. ( |A1| ) in Appendix [Aj. The larger ex- 
change energy in the DLG case, in which carriers are sep- 
arated between two different layers, has an origin similar 
to the well known increase in the exchange energy of an 
electron gas when two different spin states are occupied. 
The correlation energy, which is negative, is dramatically 
larger in DLG and strongly influenced by inter-layer in- 
teractions. These correlation contributions to the energy 
suggest that exchange-only approximations overestimate 
the degree to which screening is enhanced by beyond- 
RPA corrections. When correlation are included, the in- 
teraction energy preference for unequal partitioning of 
density between layers is also strongly reduced. This re- 
sult has an origin similar to the well known result that 
exchange-only theories badly overestimate the tendency 
of interactions to favor enhanced spin polarization in or- 
dinary electron gases. 

Finally, in Fig. [4] we present a typical result for the 
dependence of the equilibrium densities in top and bot- 



tom layers, n\ and n^, on the work function of the buffer 
layer. We clearly see that the equilibrium density in the 
top layer n\ is substantially smaller than the density in 
the bottom layer, the ratio n\/n\ changing between 1.5 
and 3.5 in the range of of <\> values explored in this fig- 
ure. Quite importantly, note that including e-e interac- 
tions at the exchange-only level severely underestimates 
the values of the equilibrium densities, especially so in 
the top layer. Including correlation effects reduces the 
energetic cost of adding electrons to the graphene lay- 
ers. In fact we find that the increase in exchange energy 
relative to the Dirac point exchange energy and the de- 
crease in interaction energy due to correlations among 
the carriers partially compensate, leading to densities in 
the graphene layers close to those implied by a Hartree 
theory which completely ignores exchange and correla- 
tion effects. This finding is surprising in comparison to 
the ordinary two-dimensional electron gas case, in which 
exchange and correlation effects always lower the energy 
cost of adding carriers and increase charge transferred 
from an electron reservoir. 
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FIG. 4: (Color online) Electro- chemical equilibrium for 
double-layer graphene on SiC in the absence of top and bot- 
tom gates (Ei = E4 = 0). All parameters are the same as 
in Figs. |2]|3l The distance between buffer and bottom lay- 
ers is equal to 3.35 A. Panel a) Equilibrium density in the 
top layer, n\ (in units of 10 12 cm -2 ), as a function of the 
chemical potential cj) (in eV) in the buffer layer. Here we 
have reported results obtained (i) neglecting intra- and inter- 
layer e-e interactions (circles), (ii) including e-e interactions 
at the exchange-only level (squares), and (hi) including both 
exchange and correlations (triangles). Panel b) Same as in 
panel a) but for the bottom layer. Note how n\ > n\. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have presented a theoretical scheme, 
based on the random phase approximation and on the 
"GoW" theory, to deal with electron-electron interac- 
tions in multi-layer graphene Fermi liquid systems con- 
sisting in which hybridization is suppressed by relative 
rotations of honeycomb lattices. We have presented nu- 
merical results for the specific case of a double layer, ex- 
plicitly demonstrating that inter-layer Coulomb interac- 
tions substantially alter the properties of the quasipar- 
ticles in this system, with respect to those of isolated 
single- layer graphene. We have also solved numerically 
the electro-chemical equilibrium problem, finding that it 
is crucial to incorporate correlation effects when going 
beyond a free-electron description of this system. When 
both effects are included, interactions play a major role in 
lowering the interaction energy of MLG systems. Because 



of strong inter-layer carrier-carrier interactions, the total 
interaction energy is relatively insensitive to charge dis- 
tribution among the layers, which is well approximated 
by a theory that includes only kinetic and electrostatic 
energies. The unusual positive exchange energy of bilayer 
graphene acts to suppress charge transfer from the buffer 
layer. Surprisingly, when both exchange and correlation 
is included corrections to the Hartree theory of charge 
transfer from the MLG buffer layer are quite small, at 
least in the N = 2 MLG system. In analysing the physics 
of larger N MLG systems, it will be essential to rec- 
ognize that layers contribute to the systems's exchange 
and correlation energy even when their carrier density 
vanishes. Because carrier-carrier interactions are weakly 
layer dependent, MLG systems might provide an inter- 
esting approximate realization of a SU(N) electron gas 
model in which 1/N expansions are quantitatively use- 
ful. This might provide an interesting attack on a limit 
in which corrections to the GoW theory can be explored 
by a rigorous semi-classical theory. 
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Appendix A: Explicit expressions for the exchange 
and RPA correlation energies of double-layer 
graphene 

For the sake of completeness, in this Appendix we re- 
port the expressions we have used to compute the ex- 
change and RPA correlation energies of DLG. In what fol- 
lows all symbols with a bar over them denote dimension- 
less quantities. We have scaled all wavevectors with the 
Fermi wavevector /cp = V^rn evaluated at the total den- 
sity n, all frequencies with the Fermi energy £ F = ^^f, 
the Lindhard response functions with the massless 
Dirac fermion density-of-states gsp / (2ttv 2 ), and the bare 
Coulomb potentials Vw with 27re 2 //cp. Here g = 4 is the 
usual spin- valley degeneracy factor. 

It is easy to prove that the exchange energy can be 
written as 



e F a e 



U+oo />+oo 
dqqVu(q)J ' 



7T 

+°° _ r+oc _ . 

dqqV 22 (q) \ dti (q^)} , (Al) 



/o Jo 
while the RPA correlation energy as 
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^ 22 (g,iO)+^2 U; fe^)J + 2 y dqqV 12 (q) J dti 5$ 12 (q,iSi)j . (A2) 
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In Eq. (A2) we have introduced the quantities (for rea- 
sons of space we will omit to write explicitly the depen- 



dence of the following expressions on the variables q and 

in) 



and 
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(A3) 



(A4) 



where 



(A5) 



L = ( 5 «ee) 2 (VllffcZ » " f? 2 Z W) , (A6) 



and 



K = ga ee (y ll xf ) +V 22 xf ] ) ■ (A7) 



The expression for $22 can be obtained from Eq. (A3) 
by interchanging 1 <H> 2. 

Fin ally, we r emark that all the quantities which in 
Eqs. (A1)-(A2) are preceded by the sign "5" indicate 



quantities defined by differences between doped and un- 



k F p=0 



> (A8) 



doped values: for example 

5xf (q, id) = xT («, in) - xf («, ifi) 

and 

= QijfaiSi) - ^ i:j (q,in)\ kFi=0 . (A9) 

Note that the dimensionless Lindhard function in the un- 
doped limit is the same for every layer and is given by 



/Cp £=0 



8 vV + fi 2 



(A10) 
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